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Nonequilibrium  dynamics  of  arbitrary-range  Ising  models  with  decoherence: 

An  exact  analytic  solution 
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The  interplay  between  interactions  and  decoherence  in  many-body  systems  is  of  fundamental  importance 
in  quantum  physics.  In  a  step  toward  understanding  this  interplay,  we  obtain  an  exact  analytic  solution  for 
the  nonequilibrium  dynamics  of  Ising  models  with  arbitrary  couplings  (and  therefore  in  arbitrary  dimension) 
and  subject  to  local  Markovian  decoherence.  Our  solution  shows  that  decoherence  significantly  degrades  the 
nonclassical  correlations  developed  during  coherent  Ising  spin  dynamics,  which  relax  much  faster  than  predicted 
by  treating  decoherence  and  interactions  separately.  We  also  show  that  the  competition  of  decoherence  and 
interactions  induces  a  transition  from  oscillatory  to  overdamped  dynamics  that  is  absent  at  the  single-particle  or 
mean-field  level.  These  calculations  are  applicable  to  ongoing  quantum  information  and  emulation  efforts  using 
a  variety  of  atomic,  molecular,  optical,  and  solid-state  systems.  In  particular,  we  apply  our  results  to  the  NIST 
Penning  trapped-ion  experiment  and  show  that  the  current  experiment  is  capable  of  producing  entanglement 
amongst  hundreds  of  quantum  spins. 

DOI:  10.1103/PhysRevA.87.042101  PACS  number(s):  03.65.Yz,  05.50.+q,  37.10.Ty,  75.10.Pq 


I.  INTRODUCTION 

Understanding  strongly  correlated  quantum  systems  in 
the  presence  of  decoherence  is  a  fundamental  challenge  in 
modern  physics.  While  decoherence  generally  tends  to  degrade 
correlations,  it  is  now  widely  appreciated  that  it  can  also  give 
rise  to  many-body  physics  not  possible  with  strictly  coherent 
dynamics  11-3],  and  can  be  used  explicitly  for  the  creation 
of  entanglement  14-6].  Regardless  of  whether  one’s  intention 
is  to  minimize  or  to  harness  decoherence,  determining  its 
effect  on  interacting  many-body  systems  is  central  to  quan¬ 
tum  simulation  17],  quantum  information  18],  and  quantum 
metrology  19].  So  far,  this  understanding  has  been  hindered 
by  the  computational  complexity  of  numerical  techniques  for 
open  systems  and  the  scarcity  of  exact  analytic  solutions. 
Exact  solutions  for  dynamics  of  interacting  quantum  systems 
in  dimensions  greater  than  one  are  rare  even  in  the  absence 
of  decoherence,  and  to  our  knowledge  no  such  solutions  have 
been  obtained  in  the  presence  of  local  decoherence. 

The  central  result  of  this  manuscript  is  an  exact  solution, 
Eqs.  (10)-(12),  for  the  time  dependence  of  all  two-spin 
correlation  functions  in  a  system  of  spins  interacting  via 
arbitrary  Ising  couplings  in  any  dimension,  and  subject  to 
local  Markovian  decoherence.  Our  solution  is  applicable 
to  a  broad  range  of  important  quantum  systems,  including 
trapped  ions  110-12],  polar  molecules  113,14],  Rydberg  atoms 
115,16],  neutral  atoms  in  optical  cavities  117,18],  optical 
lattice  clocks  119],  superconducting  qubits  120],  quantum 
dots  121],  and  nitrogen  vacancy  centers  122].  Nonequilibrium 
dynamics  in  the  presence  of  Ising  interactions  is  highly 
nonclassical,  and  can  be  exploited  in  such  systems  for  the 
generation  of  entanglement  for  use  in  quantum  metrology 
and  quantum  information  123-28]  (for  instance,  spin-squeezed 
states  123]  or  cluster  states  129]  can  be  produced).  Here 
we  apply  our  solution  to  trapped-ion  experiments  because 
(1)  the  relative  importance  of  decoherence  and  coherently 
driven  quantum  correlations  is  controllable,  (2)  the  tunable 
long-ranged  interactions  are  generically  frustrated,  making 


large-scale  numerical  simulations  impractical,  and  (3)  these 
experiments  are  the  most  developed,  with  nonequilibrium 
dynamics  already  being  explored  110].  We  will  show  that 
entanglement — in  the  metrologically  useful  form  of  spin 
squeezing — can  be  produced  amongst  hundreds  of  spins 
under  current  experimental  conditions  in  Ref.  110].  We  also 
demonstrate  that  expected  improvements  to  this  experiment 
should  allow  for  the  production  of  macroscopic  superposition 
states  (MSSs)  in  the  near  future. 

The  structure  of  the  manuscript  is  as  follows.  In  Sec.  II 
we  describe  the  model,  which  includes  coherent  dynamics 
due  to  arbitrary  Ising  interactions,  and  an  extremely  general 
form  of  Markovian  decoherence  including  both  dephasing, 
spin  relaxation,  and  spin  excitation.  In  Sec.  Ill,  we  review 
the  quantum  trajectories  approach  to  open  quantum  systems, 
which  is  central  to  our  solution  of  the  model.  Section  IV 
describes  the  application  of  quantum  trajectories  to  computing 
various  observables  relevant  to  both  relaxation  and  growth 
of  quantum  correlations.  The  idea  is  to  (a)  calculate  these 
quantities  exactly  along  a  particular  trajectory,  and  then  (b) 
carry  out  the  averaging  over  trajectories  analytically.  As  we 
will  show,  the  resulting  expressions  reveal  that  the  dynamics 
exhibits  a  transition  from  oscillatory  to  overdamped  behavior, 
an  effect  that  is  not  present  at  the  single-particle  level  or  in  a 
mean-field  treatment  of  the  interactions.  With  these  solutions 
in  hand.  Sec.  V  considers  the  applications  to  trapped-ion 
experiments  in  general,  with  a  particular  focus  on  relaxation 
and  squeezing  in  the  recent  experiment  in  Ref.  110].  We 
conclude  in  Sec.  VI  by  considering  several  open  questions 
and  promising  directions  for  future  research. 

II.  MODEL 

We  consider  far-from-equilibrium  dynamics  of  a  long- 
ranged  Ising  Hamiltonian, 

i<j 
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EIG.  1.  (Color  online)  Schematic  illustration  of  the  various 
decoherence  processes:  Ti  spin  relaxation  with  rates  ru(i,rdu  (red 
and  blue  spin,  respectively),  and  T2  dephasing  with  rate  Fei  (green 
spin).  In  atomic  systems,  one  way  this  decoherence  can  arise  is  due 
to  spontaneous  emission  from  an  excited  level  (top  panels)  [32]. 

Here  (tz  =  x,y,z)  are  Pauli  matrices,  M  is  the  total  number 
of  spins,  and  the  subscripts  (/,  j)  are  site  indices.  Our  results 
are  valid  for  arbitrary  Jtj ,  but  given  the  relevance  to  numerous 
experiments  we  will  sometimes  consider  power-law  couplings 
Jij  =  J\ri  —  where  r/  is  the  position  of  the  /th  spin 

in  lattice  units  (^  =  3  for  polar  molecules,  ^  =  6  for  Rydberg 
atoms,  and  0  <  ^  <  3  for  trapped  ions).  In  the  presence  of  local 
decoherence,  the  most  general  Markovian  dynamics  of  the 
system  reduced  density  matrix  obeys  a  master  equation  [30], 

Tip  = -  p'Hlff)-\-V{p).  (2) 

The  effective  Hamiltonian  TYeff  and  dissipator  V  have 
contributions  from  all  jump  operators  J  G  [aj 2, 

VTOT,6-;VrV8  :  1  <  7  <  W}, 

=  Vip)  =  2Y,JpJ^>  (3) 

SiWJ  allj 

where  zb  The  jump  operators  and 

give  rise  to  spontaneous  deexcitation,  spontaneous  excitation, 
and  elastic  dephasing,  respectively  (see  Fig.  1).  We  refer  to 
the  spin-changing  processes  (a^)  as  Raman  decoherence,  and 
the  spin-preserving  processes  (a^)  as  Rayleigh  decoherence. 
In  what  follows  we  assume  an  initially  pure  and  uncorrelated 
density  operator  p(0)  =  |  V^(0))(V^(0)|,  with 

|tA(0))  =  (g)^/7(a,)|CT,}.  (4) 

J 

Here  cr^  =  zb  1  are  the  eigenvalues  of  //(I)  = 

cos(0j72)^^'^3V2^  77— 1)  =  sin(0y/2)^“^'^3V2^  arbitrary 

Oj  and  (pj . 


III.  QUANTUM  TRAJECTORIES  APPROACH 

Our  approach  to  the  solution  of  Eq.  (2)  for  the  chosen  initial 
conditions  and  arbitrary  {Fud,  Tdu,  Tei}  is  based  on  the  quantum 
trajectories  method  [31],  in  which  |  V^(0))  is  time  evolved  with 
the  effective  Hamiltonian, 

and  the  dynamics  are  interspersed  with  stochastic  applications 
of  the  jump  operators.  In  Eq.  (5)  we  have  defined  Fr  = 


7 

I - 1 - 1 - 1 - 1 - 

0  t  il  f  t2  I  ts  \  ti  \  t 

Tj  (^2  —  ti)  +  (^4  —  ts)  _  ^2)  +  (t  —  14) 

FIG.  2.  (Color  online)  Series  of  Raman  flips  of  the  spin  on  site 
j  can  be  formally  accounted  for  as  a  magnetic  held  of  strength 
IJjk/Af  acting  for  a  time  In  the  notation  deflned 

below,  this  series  of  jumps  is  represented  by  the  operator  Qj  a 
d-+ (q)d-r  (^2)0-+ (t3)d-r  (^4). 


Eud  +  Edu  and  y  =  |(Eud  —  Edu).  According  to  the  standard 
prescription  [31],  a  particular  trajectory  consists  of  a  set  of 
jump  times  {^1,^2,.-.},  which  are  selected  by  equating  the  norm 
of  the  wave  function  to  a  random  number  uniformly  distributed 
between  0  and  1.  Choosing  which  jump  operator  to  apply  at 
each  time  requires  calculating  all  expectation  values  (27^  J) . 
Because  H  is  Hermitian  and  commutes  with  all  products  J, 
it  has  no  effect  on  the  selection  of  the  jumps,  which  can 
therefore  be  obtained  for  each  spin  independently  (since  the 
anti-Hermitian  part  of  TYeff  does  not  couple  different  spins). 
We  note,  however,  that  Ti  does  not  commute  with  the  recycling 
term  V{p)  defined  in  Eq.  (3).  With  the  jump  times  and  jump 
operators  in  hand,  we  define  a  string  of  uj  (time-labeled)  jump 
operators  on  site  j  as  Qj  =  JjOj)  x  •  •  •  x  The  time 

evolution  of  the  wave  function  along  a  trajectory  is  then 


\m)  =  T 


IV^(O)), 


(6) 


where  the  time-ordering  operator  T  enforces  that  the  jump 
operators  are  interspersed  in  the  time  evolution  according  to 
their  time  labels. 

The  time  ordering  of  the  Rayleigh  jumps  can  be  ignored: 
because  [d'^^,7Yeff]  =  0  and  [drj,a^]  =  ^2a^,  all  Rayleigh 
jumps  can  be  evaluated  at  ^  =  0  (their  commutation  with 
Raman  jumps  only  affects  the  overall  sign  of  the  wave 
function).  To  the  contrary,  the  Raman  jump  operators  are 
structurally  similar  to  a  transverse  field,  and  do  not  commute 
with  T^eff;  their  time  ordering  cannot  be  so  easily  accounted 
for.  However,  imagine  that  spin  j  undergoes  a  single  Raman 
jump,  created  by  applying  at  time  t.  This  jump  operator 
not  only  flips  spin  j  into  the  up  position,  but  also  removes 
all  parts  of  the  wave  function  in  which  spin  j  pointed  up 
immediately  before  time  t.  Hence  it  is  as  this  spin  pointed 
down  before  time  t,  and  up  after  time  t.  Since  spin  j  is  always 
in  an  eigenstate  of  Gj  it  is  a  spectator  to  the  Ising  dynamics, 
but  it  does  influence  the  other  spins  via  the  Ising  coupling 
Jjk',  formally  it  acts  on  spin  k  as  an  inhomogeneous  magnetic 
field  of  strength  IJjkjAf  that  pointed  down  before  t  and  up 
after  t.  Eor  a  spin  on  site  j  that  undergoes  multiple  Raman 
processes,  the  same  reasoning  allows  us  to  treat  it  as  a  field  of 
strength  that  acts  for  a  time  Xj  =  (rj^  —  where 

rJP(down)  is  the  total  amount  of  time  that  spin  j  spends  pointing 
up(down)  along  z  (see  Eig.  2).  Hence  we  are  free  to  evaluate 
all  of  the  jump  operators  at  ^  =  0  to  give  \'^)  = 
thus  ignoring  the  time  ordering  in  Eq.  (6),  so  long  as  we  evolve 
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lij/)  with  a  modified  time-evolution  operator 


U  =  exp 


-it  I  H'  +  -  iy)a^ 


(7) 


Here  r]j  =  ^  Jjk'^k  accounts  for  the  magnetic  field  of 
all  Raman-flipped  ions,  TL'  is  obtained  from  H  by  ignoring 
the  spin-spin  couplings  to  spins  that  have  undergone  Raman 
jumps,  and  y  accounts  for  the  non-Hermitian  part  of  TYgff  that 
is  not  proportional  to  the  identity  operator. 

The  expectation  value  of  an  arbitrary  operator  O  at  the 
end  of  a  particular  quantum  trajectory  is  therefore  given 
by  (O)  =  and  formally  taking  the 

average  over  all  trajectories  (denoted  by  an  overbar)  we  have 

tr[pO]  =  (O).  Along  with  Eq.  (7),  these  constitute  a  formal 
solution  for  the  dynamics  of  any  observable.  We  now  proceed 
to  derive  closed-form  expressions  for  the  transverse  spin- 
length  and  spin- spin  correlation  functions — which  have  not 
been  derived  previously  even  in  the  absence  of  decoherence. 
These  are  the  central  result  of  this  paper. 


IV.  TRANSVERSE-SPIN  LENGTH  AND 
CORRELATION  EUNCTIONS 

Relaxation  of  the  transverse-spin  length  in  an  Ising-type 
spin  model  is  a  canonical  example  of  equilibration  in  a  closed 
quantum  system  [33-35].  In  the  present  model,  such  relaxation 
occurs  due  to  a  combination  of  the  proliferation  of  quantum 
fluctuations  and  the  equilibration  with  the  environment  (de¬ 
coherence).  Our  theoretical  treatment  allows  for  both  effects 
to  be  treated  simultaneously,  and  therefore,  in  principle,  the 
disentangling  of  these  two  physically  different  (but  conse¬ 
quentially  similar)  processes.  While  such  relaxation  does  not 
directly  indicate  quantum  correlations  or  entanglement,  in 
the  absence  of  decoherence  it  nevertheless  is  entirely  due  to 
the  buildup  of  quantum  correlations — at  the  mean-field  level 
coherent  relaxation  is  absent.  The  correlations  that  develop 
during  the  dynamics  can  be  understood  in  more  detail  by 
looking  at  two-spin  correlation  functions — for  instance,  these 
furnish  a  complete  description  of  spin  squeezing  [23].  In  the 
absence  of  decoherence  and  for  all-to-all  coupling  (^  =  0), 
it  is  well  known  [23]  that  the  transverse- spin  component 
revives  at  a  time  =  Nnh/ill),  and  that  a  highly  entangled 
macroscopic-superposition  state  (MSS)  appears  at  time  Xyjl. 
This  state  is  characterized  by  vanishing  spin  length  but 
maximum  transverse- spin  fluctuations,  and  our  solution  can 
be  used  to  assess  the  robustness  of  such  fluctuations  against 
decoherence. 

To  calculate  the  transverse- spin  length  along  a  particular 
trajectory,  we  assign  the  discrete- valued  variables  IZj  and 
Tj  to  each  lattice  site,  which  count  the  number  of  Raman 
jumps  and  Rayleigh  jumps,  respectively.  As  we  have  dis¬ 
cussed  earlier,  all  jump  operators  can  be  applied  at  ^  =  0, 
and  therefore  specifying  the  [1Zj,Tj,Tj}  on  each  site  fully 
determines  transverse- spin  length  along  that  trajectory.  The 
transverse  spin  component  in  direction  (p  and  on  site  j 
is  given  simply  in  terms  of  the  spin-raising  operator  on 
that  site,  (5'J)  =  cos(^(5'y)  -h  sin(^(5'J)  =  Re[^“^‘^(d'^^)],  and 
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generalizing  solutions  obtained  in  Refs.  [33,36],  we  find 

+  _  otjPj  gk[2akt{y  -  iJjk/Af)] 

2gji2yt)  gkilytuk) 

(8) 

Here  gy(x)  =  aj=Snjfi  being  the 

Kronecker-delta  symbol),  fij  =  and  the  details  of  the 

calculation  are  given  in  Appendix  A.  Defining  a  function 
that  determines  the  probability  distribution  of  these 
variables  on  a  given  lattice  site,  we  have 

=  I  ■■■  I  Y\dTkmik,J^k,r,){a+).  (9) 

all7^  ^  ^  k 

Equation  (9)  constitutes  a  formal  solution  for  and  it 
can  always  be  evaluated  efficiently  by  averaging  {drj~)  over 
stochastically  generated  trajectories.  However,  because  the 
noise  is  uncorrelated  from  site  to  site,  and  accordingly  the 
expression  inside  the  product  of  Eq.  (8)  depends  only  on  local 
stochastic  variables,  these  sums  and  integrals  factor  into  AT 
independent  sets  (each  over  the  three  stochastic  variables), 
admitting  closed-form  expressions.  We  note  that,  while  most 
quantum  trajectories  simulations  average  over  a  finite  number 
of  stochastically  selected  trajectories,  thus  giving  rise  to 
statistical  errors,  our  approach  effectively  averages  over  an 
infinite  number  of  trajectories.  Hence  our  solutions  have  no 
statistical  approximation.  We  have  checked  for  small  spin 
systems  that  they  produce  results  identical  to  a  full  numerical 
solution  of  the  master  equation  [Eq.  (2)]. 

At  this  point,  to  avoid  unnecessary  complications  in  the 
ensuing  expressions,  we  will  take  our  initial  state  to  point  along 
the  V  axis  (0  =  7r/2,  ^  =  0),  but  our  results  easily  generalize 

[37] .  In  Appendix  B  we  explain  how  to  evaluate  these  sums 
and  integrals,  and  here  we  simply  quote  the  result.  Defining 

O(/,0  =  e~^^[cos(ty/s^  —  r)  -h  kt  sinc(^V5^2^^)], 

with  k  =  Vy/l,  s  =  2iy  2J/M,  and  r  =  EudEdu,  we  find 

n  (10) 

k^j 

where  the  total  decoherence  rate  is  defined  E  =  f  (Er  -h  Egi) 

[38] .  Similar  calculations  to  those  described  above  yield  spin- 
spin  correlation  functions 

n  (11) 

U  (12) 

itUM 

with 

=  e~^\is  —  2y)t  sinc(^\/ —  r)  (13) 

and  11,  V  =  zb.  These  correlation  functions,  along  with  similar 
ones  obtained  by  interchange  of  the  site  indices,  completely 
determine  the  spin- spin  correlations.  Each  instance  of  or 

in  an  observable  generates  an  overall  multiplicative  factor  of 
e~^\  which  would  also  occur  if  the  decoherence  were  treated 
in  the  absence  of  interactions;  this  is  the  effect  of  decoherence 
at  the  single-particle  level.  The  structure  of  0(/,Q  captures 
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EIG.  3.  (Color  online)  Plots  of  O(/,0  for  y  =  0  and  Fr/F^  e 
{0, 1/4, 1/2, 3/4, 1},  showing  a  transition  from  oscillatory  to  damped 
behavior. 

the  interplay  of  decoherence  and  the  many-body  physics,  and 
could  not  have  been  deduced  without  our  exact  treatment.  The 
additional  damping  present  in  O(/,0  has  a  straightforward 
physical  origin,  which  is  made  particularly  clear  by  the 
quantum  trajectories  formalism.  Even  if  a  spin  on  site  j  has 
not  undergone  any  quantum  jumps,  it  experiences  a  fluctuating 
(in  time  and  space)  longitudinal  magnetic  field  r]j  due  to  the 
various  other  spins  that  have  undergone  quantum  jumps,  and 
thus  experiences  additional  dephasing. 

We  note  that  for  y  =  0  the  function  undergoes 

a  qualitative  transition  from  oscillatory  (v^  >  r)  to  damped 
(s^  <  r)  behavior  when  Vf  =  =  4J /J^  (Fig.  3).  Therefore, 

when  there  is  only  one  coupling  strength  /,  for  instance  in  the 
case  of  all-to-all  or  nearest-neighbor  couplings,  dynamics  of 
the  transverse- spin  length  and  correlation  functions  undergoes 
the  same  transition.  We  note  that  the  factors  of  in 
Eqs.  (10)-(12)  may  make  this  transition  difficult  to  observe 
experimentally,  since  the  correlations  are  rapidly  suppressed 
at  the  critical  F^.  It  is  important  to  realize  that,  while  such 
a  transition  is  expected  for  a  single  spin  in  the  presence  of  a 
coherence  restoring  drive  (i.e.,  a  transverse  field),  it  does  not 
occur  for  a  noninteracting  system  in  the  absence  of  a  transverse 
field.  Hence,  in  the  Ising  model,  this  transition  is  a  manifestly 
many-body  effect,  arising  from  the  competition  of  decoherence 
and  the  coherent  dynamics  generated  by  the  interactions.  It 
is  interesting  to  note  that,  if  the  Ising  interaction  had  been 
treated  at  the  mean-field  level,  the  system  would  behave  as 
a  collection  of  independent  spins  undergoing  decoherence 
and  experiencing  a  (self-consistently  determined)  longitudinal 
field.  Thus  this  transition  cannot  be  captured  at  the  mean-field 
level. 


V.  APPLICATION  TO  A  TRAPPED-ION 
QUANTUM  SIMULATOR 

Trapped-ion  systems  can  simulate  the  Hamiltonian  in 
Eq.  (1),  and  can  accurately  measure  the  decoherence  rates 
Fei,  Fud,  and  F^u-  We  note  that  both  the  Bom-Markov  ap¬ 
proximation  and  the  assumption  of  uncorrelated  decoherence 
processes  are  extremely  well  justified  for  trapped-ion  systems 
[32].  Sample  averaged  spin-length  and  spin-spin  correlation 
functions  are  easily  measured  in  these  experiments  by  looking 
at  the  length  (and  its  shot-to-shot  fluctuations)  of  various 
projections  of  the  Bloch  vector.  In  the  trapped-ion  experiments 
discussed  in  Ref.  [10],  when  ^  =  0  the  time  scale  at  which 


FIG.  4.  (Color  online)  (a)  Example  of  how  spin  squeezing  is 
affected  by  decoherence  (dashed  lines  are  F  =  0.06/;  solid  lines 
are  for  F  =  0)  for  long-ranged  (red,  ^  =  0)  and  short-ranged  (blue, 
^  =  3)  interactions. 


quantum  correlations  become  important  for  these  observables, 
Tc,  scales  with  some  power  of  the  ion  number:  ~ 

for  spin  squeezing,  for  transverse- spin  relaxation,  and 

for  the  creation  of  MSS’s.  Taking  A/'  =  100  and  F  = 
0.06/,  as  is  typical  in  that  experiment,  we  expect  the  proper 
incorporation  of  decoherence  to  be  quantitatively  important 
even  for  spin  squeezing,  despite  it  being  a  relatively  short-time 
indication  of  entanglement. 

Equations  (11)  and  (12)  allows  us  to  exactly  calculate  the 
effect  of  decoherence  and  the  finite  range  of  interactions  on  the 
maximum  spin  squeezing  achievable  in  experiment.  Eigure  4 
shows  the  expected  maximal  squeezing  and  antisqueezing  as 
a  function  of  angle  xj/  in  the  yz  plane  for  ^  =0,3,  with  F  = 
0.06/  and  Fgi  =  8Fud  =  8Fdu  (typical  experimental  numbers 
in  [10]).  The  effects  of  decoherence  are  more  pronounced  for 
shorter-range  interactions  due  to  the  longer  time  scales  for 
maximal  squeezing.  For  this  calculation  the  spins  are  assumed 
to  be  initialized  (prior  to  the  Ising  dynamics)  in  a  coherent  state 
pointing  along  the  x  axis,  and  we  define  the  spin- squeezing 
parameter. 


VJ/AS'l' 

W) 


(14) 


with  AA  =  S'^  =  ^  V’  V  “ 

cos(\j/)aJ  -h  sin(\l/)drj.  Notably,  in  the  long-ranged  interaction 
limit  [39]  the  spin  squeezing  is  degraded  but  still  quite 
appreciable  under  current  experimental  conditions  (for  short- 
range  interactions,  ^  =  3,  the  spin  squeezing  is  completely 
destroyed  by  decoherence). 

For  the  current  experimental  parameters  there  is  essen¬ 
tially  no  spin  revival,  and  no  indication  of  a  MSS  at  rr/2. 
However,  assuming  97  ions  and  expected  improvements  in 
the  experiment  [10]  (a  roughly  50-fold  increase  in  the  ratio 
//  F),  in  Fig.  5  we  show  that  transverse- spin  revivals  begin  to 
appear.  In  that  figure,  the  dotted  line  is  obtained  by  treating 
the  decoherence  at  the  single-particle  level,  which  amounts 
to  attaching  a  decaying  exponential  to  the  operators 
and  .  We  note  that  this  treatment  would  be  exact  if  the 
decoherence  were  only  of  the  Rayleigh  type  (Fj-  =  0).  The 
solid  line  is  the  full  solution  from  Eq.  (10);  the  large  (~35-fold) 
discrepancy  between  the  single-particle  and  exact  results 
indicates  that  a  proper  accounting  of  Raman  decoherence 
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FIG.  5.  (Color  online)  Transverse- spin  relaxation  and  revivals 
for  ^  =  0,  with  parameters  corresponding  to  expected  experimental 
capabilities  in  [10]  (blue  solid  line).  The  dotted  blue  line  is  obtained  by 
treating  decoherence  at  the  single-particle  level,  and  underestimates 
the  detrimental  effect  of  Raman  decoherence  by  about  a  factor  of  35. 
Inset:  transverse- spin  fluctuations  peaking  at  time  Tr/2.  Experimental 
parameters,  exact  treatment  (red  solid  line);  experimental  parameters, 
single-particle  treatment  (red-dotted  line);  no  decoherence  (black 
dashed  line). 

processes,  which  is  achieved  in  our  solution,  will  be  essential 
for  understanding  the  behavior  of  this  experiment.  In  the  inset 
of  Fig.  5  we  plot  the  transverse- spin  fluctuations  AS^  at  times 
near  Ty/2.  The  peak  in  these  fluctuations  (which  would  achieve 
unity  in  the  absence  of  decoherence)  is  a  result  of  strong 
transverse-spin  correlations  in  an  emerging  MSS,  indicating 
that  the  expected  improvements  to  the  experiments  will  bring 
within  reach  the  production  of  MSS’ s  of  "^100  ions  in  the  near 
future. 


VI.  CONCLUSIONS 

These  calculations  provide  a  rare  glimpse  into  the  exact 
structure  of  relaxation  dynamics  in  an  open  and  strongly 
interacting  quantum  system.  The  surprisingly  large  detrimen¬ 
tal  effect  of  Raman  decoherence — which  far  exceeds  that 
expected  from  a  single-particle  picture — is  the  result  of  an 
interaction-mediated  back  action  of  decohered  spins  on  the 
remaining  many-body  system,  in  which  each  Raman-flipped 
spin  behaves  as  a  (temporally)  fluctuating  magnetic  field. 
Interestingly,  our  exact  solution  reveals  that  the  net  effect 
of  this  back  action  is  (for  ^  =  0)  a  global  rotation  in  spin 
space,  which  in  principle  can  be  corrected  if  the  arrival  times 
of  spontaneously  scattered  photons  are  recorded.  These  ideas 
will  be  explored  in  more  detail  in  future  work. 


Note  added.  Recently,  we  were  informed  of  simultaneous 
calculation  of  spin- spin  correlation  functions  for  quantum 
Ising  models  in  the  absence  of  decoherence  [40] . 
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APPENDIX  A:  EXPRESSIONS  EOR  SPIN  LENGTH 
AND  CORRELATION  FUNCTIONS  ALONG 
A  SINGLE  TRAJECTORY 

As  in  the  main  text,  we  write  the  state  of  a  single  spin 
as  where  is  an  index  that  takes  on 

the  eigenvalues  of  the  operator  fj{\)  =  co^{0jl2)e^^}/'^, 
and  //(—I)  =  ^m{0jl2)e~^'^j/^.  The  initial  state  of  the  entire 
system  is  taken  to  be  a  direct  product  of  states  for  each 
individual  spin: 

1^(0))  ^ 

j 

=  E  /iK)  X  •••  X 


As  discussed  in  the  main  text,  evaluating  all  of  the  Rayleigh 
jumps  at  ^  =  0  can  be  accomplished  by  changing  cpj  cpj  -h 
7t!Fj,  which  rotates  the  spin  on  site  j  by  an  angle  tt  if  it 
has  undergone  an  odd  number  of  Rayleigh  jumps.  Raman 
jumps  can  be  incorporated  at  ^  =  0  by  setting  Oj  =  0(7r)  if 
the  final  Raman  jump  on  site  j  was  We  therefore 

define  fjicf])  to  be  the  modification  of  under  those 

transformations.  We  also  note  that  spins  having  undergone 
one  or  more  Raman  jumps  are  treated  as  an  effective  magnetic 
field,  and  not  included  in  the  spin- spin  coupling  term  of  the 
Hamiltonian.  This  is  accomplished  by  changing  Jfj  ataj  Jtj 
(with  otj  =  Suj  ,o)  in  the  Hamiltonian  and  including  an  effective 
magnetic  field.  Therefore,  we  can  write  the  (time-dependent) 
wave  function  evolving  under  U  [as  defined  in  Eq.  (6) 

of  the  manuscript]  as 


j 


/i(^0 


•  X 


or, . . . 


Nh 


(A2) 


In  order  to  calculate  the  transverse- spin  length  we  will  first  calculate  {aj~)  =  \  ,  and  at  the  end  obtain 

(5")=Re(E,-(^;». 

Let’s  imagine,  in  particular,  calculating  1^(0)  (there  is  nothing  special  about  the  first  spin;  this  just  makes  the  notation 

in  what  follows  less  confusing).  Because  the  wave  function  enters  twice,  this  would  involve  two  sums  like  the  one  in  Eq.  (A2), 
over  Gj  and  Gj^,  but  very  few  terms  survive:  we  need  g^  =  —I,  =  I,  and  for  all  j  I  we  must  have  Gj  =  Gj^,  so  the  matrix 
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element  is  given  by 


=  /r(l)/i(-l)  Y!  |/2(o^2)P  X  •••  X  |Xv'(CT_^)|^exp 

If  aj  =  0,  the  jth  spin  always  has  a  well-defined  value  of  and  the  choice  to  include  the  term  ya^j  in  the  exponentiated 
sum  (or  not)  only  affects  the  overall  normalization  of  the  wave  function.  By  multiplying  y  by  aj  (in  the  exponentiated  sum), 
we  have  chosen  to  not  include  the  term  ycfp  and  this  is  properly  accounted  for  when  normalizing  this  expectation  value 
below.  In  order  to  obtain  we  must  divide  by  the  (nonconserved)  normalization  of  the  wave  function  .  Defining 

gj{x)  =  (as  in  the  manuscript),  we  obtain 


(A3) 


(^r)  = 


/r(i)/i(-i) 

^i(2k0 


E 


g2{2a2Yt) 


gM(2aMyt) 


exp 


2it 


/  V  1 


(A4) 


This  expression  can  be  simplified  by  making  the  following  set  of  observations:  (1)  /*(l)/^(— 1)  =  sin(0i)/2  [as  in  the 

text,  pj  =  (—1)*^^  ],  (2)  the  a\  in  the  exponent  is  irrelevant,  because  if  it  takes  the  value  0  the  entire  expression  for(i^(0|dj^  1^(0) 
vanishes,  and  (3)  the  summand  factorizes  into  a  product  where  each  term  contains  only  local  (i.e.,  on  a  single  site)  variables,  and 
hence  the  sum  of  products  can  be  exchanged  for  a  product  of  sums.  Taking  observations  (l)-(3)  into  account,  we  obtain 


(^1+) 


sin(0i)^^‘^i 

sin(0i)e^‘^i 

2gi(2Kf) 


z\\2 


“1^1  n  E 


i=2 

u 


gji2ajyt) 


exp[2itaj(Tj(Jij /Af  + 


(XiPi  H 


j=2 


.2iJ,<z,/M8j\2ajt(y  -  iJij/Af)] 
gj(2ytaj) 


Therefore,  we  can  write 


M 

(5">=ReY 

y=i 


sin(0y)^^‘^i 

2gj(2yt) 


n  (' 

k^j  ^ 


.2iJ:,zUMSk{2Mkt{y  -  iJjk/X')] 

gk{2ytak) 


(A5a) 

(A5b) 


(A6) 


The  calculation  of  correlation  functions  follows  from  extremely  similar  considerations.  For  instance,  let’s  consider  ^Jk^ 
{oTj  OT]^).  In  this  case,  the  operators  in  the  expectation  value  only  connect  two  states  if  —a'j^  =  aj  =  I,  =  I,  and 

cr/^  =  whenever  I  ^  [j,k}.  Therefore,  much  as  before  we  have 

-  Agj{2yt)gki2yt)  ^  gi(2yta,) 


Computing  correlation  functions  involving  a  single  d^,  such  as  =  (d^^d^),  can  be  achieved  by  inserting  into  the  sum  in 
Eq.  (A4),  yielding 


^m{0k)e^^'^akPk 


Qfk  = 


2gk(2yt) 


cos^(0/2)e 


—2yt  _ 


sm 


'(e/2)e^y‘ 


gj{2yt) 


+  ( 1  —  a  j)K 


rr  (^2i Jt, z,/MSi\.2ait(y  -  iJki/Af)]\ 
^  ^  '  gi(2ytai)  ) 


J  1  1  \  ' 

^  HUM  ^ 


(A8) 


Assuming  one  or  more  Raman  flip  occurred,  the  variable  kj  takes  on  the  values  ±1  if  the  final  Raman  jump  is 


APPENDIX  B:  ANALYTIC  EVALUATION  OE  STOCHASTIC 
AVERAGING  OE  TRAJECTORIES 

At  this  point  in  the  calculation,  for  clarity  of  presentation, 
we  set  (pj  =  0  and  Oj  =  7t  12  (for  all  j),  so  all  spins  point  along 
the  V  axis  at  ^  =  0.  Defining  V(1Z,!F,t)  to  be  the  probability 
distribution  of  the  variables  IZ,  T,  and  r  on  a  single  site  (it  is 


the  same  on  every  site),  the  trajectory  averaged  expectation 
value  is  given  by 

^=EE/  dxi  . . .  f 

all7^  allJ^  J  J  y. 

(Bl) 
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To  begin,  we  note  that  the  probability  distribution  can  be 
decomposed  as  which  is  valid 

because  the  probability  of  Rayleigh  jump  is  independent 
of  whether  a  Raman  jump  has  occurred  (and  vice  versa). 
The  occurrence  of  random  processes  follows  a  Poissonian 
distribution,  so  Vq\(J^)  =  /4)^ /T\,  and  we  have 

also  calculated  Vr.  The  result  depends  on  whether  IZ  is 
even  or  odd  (the  proof  is  simple  but  requires  some  careful 
reasoning,  and  we  do  not  give  it  here).  We  parametrize 
the  7^-odd  solution  by/x  =  (7^  —  l)/2  (which  will  run  over 
all  non-negative  integers),  and  we  parametrize  the  7^-even 
solutions  by  /X  =  (7^  —  2)/2  (once  again  running  /x  over  all 


non-negative  integers),  and  obtain 


By  evaluating  the  sum 


,-rel(/2 


(B2) 


(B3) 


JF=0 


we  obtain 


J 


(^n  = 


^-reiU2 

2  cosh(2}/0 

g-reiU2 

2  cosh(2}/0 


E  J 


dxM 


jrJ  coshdrla,) 


It,  ■'  cosh(2Kt«,) 

Because  a j  gives  1  if  there  have  not  been  any  Raman  flips  at  site  j  and  0  otherwise,  the  term  in  the  first  square  bracket  is  just  the 
probability  that  there  has  been  no  Raman  flip  on  site  j,  which  is  given  by  ^“^rV2  cosh(2]/Q  [this  comes  from  evolving  the  wave 
function  of  a  single  spin  pointing  along  x  with  the  effective  Hamiltonian  Eq.  (4)].  Therefore,  we  have 


1 

2^ 


-r« 


j^k  TZk 


(B4) 


Defining  s  =  2iy  -\-  2J /M,  the  trick  is  now  to  evaluate  the 
quantity 


OO  . 

_  ^-rrV2  cosh(/5'0  +  ^  /  dr  . 

7^=l 


(B5) 


The  second  equality  follows  from  pulling  off  the  7^  =  0  term  in 
the  sum,  which  represents  the  probability  of  having  no  Raman 
flip  (and  so  we  can  set  r  =  0  and  a  =  1  in  this  term).  The 
second  term  represents  the  probability  for  any  finite  number 
of  Raman  flips,  and  hence  we  must  keep  r  arbitrary  but  can  set 
O'  =  0.  The  integral  over  r  can  be  evaluated  by  using  the 
identity 

f  dx{f  -  =  (2f)^+*  jidixt)ii\  ^ 

J-t  (xT 

where  j  is  a  spherical  Bessel  function.  The  remaining  sum 
over  IZ  can  be  recognized  as  a  generating  function  for  the 


spherical  Bessel  functions  (or  a  derivative  thereof).  Defining 
parameters  k  =  r^/2  and  r  =  Vud^du,  and  functions 


F(x,y)  =  sinc(yx^  —  y), 

cos(y/x'^  —  y)  —  cos(x) 


G(x,y)  = 


we  obtain 


(B7) 

(B8) 


c^(J,t)  =  e  cos(st) kt  e  F(st,rt) st  e  ^^G(st,rt) 
=  e~^^[cos(ty/s^  —  r)  -h  kt  sinc(ty/s^  —  r)]. 

(B9) 

We  can  now  write  out  the  exact  solution, 

W)  =  ^Re  n  (BIO) 

j  k^j 

Because  Eqs.  (A6)  and  (A7)  have  such  a  similar  structure, 
the  stochastic  averaging  of  correlation  functions  is  almost 
identical,  and  leads  to  the  similar  expressions  given  in  the 
main  text  [Eqs.  (10)  and  (11)]. 
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